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IX  AISTRACT  (Mtjomwn  JOOyovmi 

The  design  of  a  system  or  circuit  in  which  stability  is  a  key  constraint  frequently  leads  to  an  optimization 
problem  over  the  space  of  functions  analytic  on  the  right  half  plane  (R.H.P.)  Mathematical  techniques  for  solving 
such  optimization  problems  for  mean  square  error  (L*  error)  criteria  have  been  widespread  in  engineering  since 
the  time  of  Wiener.  Much  of  this  research  goes  to  developing  techniques  for  handling  worst  case  error  {L°°  error) 
criteria.  These  occur  naturally  in  design  of  control  s^'stems  and  amplifiers.  Practically  speaking  there  is  evidence 
that  frequency  domain  L°°  criteria  control  system  designs  have  desirable  robustness  properties.  The  ultimate 
objective  is  to  develop  a  new  CAD  approach  to  MIMO  control  design  which  has  the  flavor  of  classical  control  as 
well  as  a  systematic  approach  to  worst  case  frequency  domain  design  as  it  occurs  in  many  areas.  The  promise  of 
this  approach  is  sufficient  to  have  attracted  many  investigators  and  it  is  currently  the  focus  of  much  attention.  « 

This  research  addresses  many  aspects  of  the  problem.  They  range  from  the  development  of  computer  al¬ 
gorithms  of  a  radically  different  type  to  the  discovery  of  theoretical  methods  for  understanding  computational 
design.  Also  considerable  progress  was  made  in  extending  existing  H°°  control  to  nonlinear  plants. 

Another  major  effort  involves  computer  algebra  for  systems  research.  The  objective  is  to  treat  (on  a  computer) 
systems  formulas  of  the  type  an  investigator  would  manipulate  by  hand.  Considerable  software  was  develop>ed 
along  these  lines. 
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FINAL  REPORT 

The  research  concerns  control  but  focuses  on  substantially  different  parts  of 
the  subject,  namely,  nonlinear  systems,  optimization  theory  and  algorithms  for  fre¬ 
quency  domain  design,  and  computer  algebra  tailored  to  systems  and  control  research. 


Nonlinear  systems 

The  modern  approach  to  worst  case  design  in  the  frequency  domain  arose  from 
studies  of  amplifier  design  the  "dual”  problem  of  making  a  circuit  dissipative  using 
feedback.  For  linear  systems  key  cases  of  this  were  solved  in  1965  (SISO)  by  Youla  and 
Saito  and  (MIMO)  in  1976  by  Helton.  In  the  early  80’s  Zames  and  Francis  formulated 
H°°  control  and  solved  the  math  problem  by  drawing  on  the  earlier  solutions  to  this 
circuits  problem.  In  the  beginning  the  subject  of  control  evolved  quickly  in 
significant  part  because  key  math  problems  were  already  reasonably  understood  by 
operator  theorists.  I  participated  in  this  earlier  work  (e.g.  solved  the  MIMO 
control  problem  with  Zames  and  Francis,  also  Pearson  and  Chang)  but  at  the  same 
time  begin  pushing  in  new  directions:  nonlinear  plants  and  an  H’**  approach  to 
classical  control. 
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Of  the  various  solutions  to  CTRL  one  which  is  easy  to  implement  and  numeri¬ 
cally  sound  is  the  Doyle-Glover-Kargonekar-FVands  DGKF  two  Riccati  equation  solu¬ 
tion.  Consequently  extending  this  to  nonlinear  plants  is  of  considerable  importance. 
Through  the  last  3  years  there  has  been  considerable  progress  by  Isidori  and  coworkers 
and  by  our  group  (Ball  Helton  Walker  Zhan).  Isidori  et  al  find  local  sufficient  condi¬ 
tions  and  compute  (with  Krener’s  software)  power  series  solutions  to  model  problems. 
All  of  these  approaches  assume  something  like  the  dimension  of  the  compensator’s 
state-space  equals  that  of  the  controller  state-space. 

Our  results  [BHW]  say: 

Result  1  The  DGKF  equations  are  2  equations  each  in  n  dimensional  space.  These 
generalize  to  the  nonlinear  case  as  one  Riccati  P.D.E.  in  2n  variables.  A  positive 
solution  to  it  is  necessary  for  CTRL  to  have  a  solution.  In  the  linear  case  the  2n 
dimensional  Riccati  P.D.E.  easily  yields  the  two  DGKF  equations. 

There  is  a  surprisingly  strong  yet  general  separation  theorem  which  limits  the  type 
of  controller  you  can  effectively  use. 

There  is  a  formula  which  is  reasonable  to  try  for  the  controller.  Only  the  input 
term  for  the  controller  is  a  compromise. 

Wei  Zhan  and  I  analyzed  completely  those  compensators  whose  state-space  have 
the  same  dimension  as  the  plant’s  state-space: 

Result  2  There  is  a  strange  type  of  equation  which  is  a  mixture  of  a  first  order 
P.D.E.  and  best  approximation  operators  which  we  call  a  Tchebychev  Riccati  PDE. 
Existence  of  a  positive  solution  to  this  equation  is  equivalent  to  a  type  of  solution  to 
CTRL. 


It  is  extremely  unlikely  that  Tchebychev  Riccati  PDE  will  ever  be  solved  exactly. 
However,  now  we  know  the  enemy  and  this  should  help  organize  compromises  in  a 
systematic  fashion. 

Also  there  is  progress  on  evaluating  performance  of  piecewise  linear  systems.  We 
took  a  typical  architecture  (a  la  Campo  Morari)  for  a  system  with  saturation  and 
extracted  one  of  the  key  computational  difficulties.  These  systems  are  piecewise  ^ — 
linear  and  continuous.  Work  in  progress  with  Ball  shows  that  a  key  object  for  a 
dissipative  system  ,  called  a  storage  function,  must  be  continuous.  We  then  made  a 
natural  compromise.  The  continuity  of  the  storage  function  forces  constraints  which  .ion_ 
make  analyzing  such  systems  not  a  Linear  Matrix  Inequality.  We  found  a  sequence - 
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of  steps  which  extracted  the  non  LMI  part  and  allowed  one  to  solve  the  problem  of 
determining  performance  of  such  systems  by  doing  first  an  LMI  check,  then  a  side 
test  then  an  LMI,  etc. 

There  are  still  basic  theoretical  issues.  In  the  last  year  James  and  Baras  have 
necessary  and  sufficient  conditions  on  the  control  problem  when  one  allows  an 
infinite  dimensional  state-space.  Under  a  saddle  point  assumption  these  reduce  to 
one  due  to  van  der  Schaft.  Krener  has  results  of  a  similar  tone.  James  will  visit  here 
for  most  of  the  spring  quarter. 

A  main  open  issue  in  the  subject  is  finding  compromise  solutions  to  the  equations 
produced  above.  Ultimately  I  see  much  of  the  field  as  consisting  of  sensible  ways  of 
finding  conservative  solutions  to  the  equations  which  arise  in  the  theoretical  studies 
above. 

Optimization  over 

Much  of  my  effort  goes  to  studying  a  basic  question  of  worst  case  frequency  domain 
design  where  stability  of  the  system  is  the  key  constraint.  This  is  the  optimization 

problem  which  is  crucial  in  several  branches  of  engineering. 

The  fundamental  H°°  problem  of  control. 

First  we  state  the  core  mathematics  problem  graphically.  At  each  frequency  ut  we 
are  given  a  set  5u,(c)  C  C^,  called  the  specification  set  The  objective  is  to  find  a 
function  T  with  no  poles  in  the  R.H.P.  so  that  each  T{jui)  belongs  to  5w(c).  In  fact 
there  is  a  simple  picture  to  think  of  in  connection  with  a  design 


Typically  there  is  a  nested  family  of  tarf^t  sets  5w(c)  parameterized  by  a  perfor¬ 
mance  level  c;  the  smaller  the  sets  the  better  the  performance.  For  the  optimal  c  a 
solution  T  exists  but  no  solution  exists  for  tighter  specs. 

The  Horowitz  templates  of  control  can  be  transformed  into  this  type  of  picture. 
When  each  ^^(c)  is  a  ‘‘disk”  this  problem  is  solved  by  transformations  of  “classical 


3 


pure”  mathematics  done  in  the  late  1970*8  by  Helton.  Many  different  solutions  to 
this  problem  in  many  different  coordinates  were  worked  out  by  engineers  in  the  last 
15  years  since  it  is  the  subject  of  control.  Competing  constraints  and  plant 
uncertainty  lead  inunediately  to  spec  sets  which  are  not  disks. 

The  graphical  problem  of  Fig2  can  be  formulated  analytically  in  terms  of  a  per¬ 
formance  function  F  as 

(OPT)  Given  a  positive  valued  function  F  on  R  x  (which  is  a  performance  mea¬ 
sure),  find  7*  >  0  and  /*  in  As  which  solve 

7*  =  in/  snp  r(w,/(iu;)) . 


and  this  of  course  is  what  one  puts  in  a  computer.  Collaborators  and  I  have  a  very 
broad  based  attack  on  the  problem  which  addresses  most  aspects  of  it. 

fVom  qualitative  theory  to  numerical  algorithms  and  diagnostics. 

While  little  was  known  about  this  problem  10  years  ago  there  has  been  a  lot  of 
progress,  and  now  we  have  substantial  theory.  In  progress  is  an  elementary  book 
with  Merino  on  control  system  design  which  gives  our  methodology  for  setting  design 
problems  as  formal  optimization  problems.  Then  our  software  solves  the  problems. 
The  software  runs  under  Mathematica  and  can  be  obtained  from  anoptQucsd.edu. 

We  shall  not  sketch  all  that  is  known  about  OPT  but  mention  one  dramatic 
qualitative  result  (with  D.  Marshall  ) 

Result  3  For  a  **properly  formulated”  SISO  control  problem  the  optimal  compensator 
is  unique. 

Here  no  convexity  is  assumed. 

Ironically  one  of  the  most  practical  results  on  an  optimization  problem  is  charac¬ 
terization  of  the  optimum,  since  this  is  the  basis  for  numerics.  Our  result  is  easier 
to  state  on  the  unit  disk  A  and  the  unit  circle  T  rather  than  on  the  R.H.P.  and 
the  ju-axis.  Also  we  state  it  only  for  the  N=:2  MIMO  case.  Roughly  the  optimality 
condition  for  solutions  to  OPT  is 

Result  4  Given  F  a  smooth  function.  Necessary  an  sufficient  for  a  smooth  function 
T*  in  H'^  satisfying  a(e**)  =  (c*^,  r*(e^))  is  never  0  on  T  to  be  a  local  solution  to 

OPT  is 


4 


/  r(e**,r*(e’*))  is  constant  in  e**. 

II  There  exist  Fi  and  Fj  analytic  on  the  disk  and  A  a  positive  function  on  the 
circle  such  that  for  all  €  T, 

^  («",r(«"))  =  e»A(e‘*)F,(e’') 

III  A  condition  on  second  derivatives  ofT. 

A  typical  computational  strategy  is  to  apply  Newton’s  method  to  (I)  and  (II)  above 
to  solve  them  thereby  solving  OPT.  Even  when  the  spec  sets  were  disks,  subcases  of 
which  have  been  studied  for  80  years,  Newton’s  method  was  never  succcssfLlly  applied 
directly  to  this  problem.  The  difficulty  is  that  the  problem  is  highly  degenerate. 
However,  recently  O.  Merino,  T.  Walker  and  I  obtiuned 

Result  S  There  is  a  functional  analysis  transform  of  (I)  and  (II)  which  yields  nonde¬ 
generate  equations  and  so  Newton's  method  applied  to  them  is  second  order  convergent. 

Consequently  we  are  finally  obtaining  satisfactory  computational  methods  for  solv¬ 
ing  OPT.  In  order  to  show  this  we  give  an  example  and  optimize  it  using  Newton 
iteration  (Table  1)  and  the  Disk  iteration  (Table  2). 

EXAMPLE  The  problem  is  to  solve  OPT  for  the  performance  function 

r(c’^,  Zi,  Zj)  =  |zi|*+|z3|*-|-|100+c'*Zi-l-.l(Zi-|-Za+ZiZ3)|*-|-|100-|-e’*Z3-|-.l(Zi-l-Z2+ZiZ2)|* 

The  (absolute)  optimal  value  is  7*  =  3800,  attained  at  two  different  points  in  function 
space:  the  constant  functions  f^  —  (30,  —30)  and  f^  =  (—30,30).  There  are  no  other 
local  solutions  to  OPT. 

The  meaning  of  each  column  in  Table  1  is  as  follows:  It  (Iteration  number).  Value 
(Current  value  of  sup|r(*,/)),  ||  f  -  f*  ||  (The  true  error,  f*  is  the  solution)  ,  OTl 
(checks  for  equation  I  oi  Result  4),  OT2  (^ecks  for  equation  II  of  Result  4),  NED  ( 

A  measure  of  numerical  error). 

The  discretization  of  the  problem  is  carried  out  by  sampling  fimctions  on  a  grid 
of  256  equally  spaced  pmnts  on  the  unit  circle.  The  Newton  iteration  is  initialized  at 
/o(e^)  =  (29.6  -f  .le^,  —30.4  —  .OOOle^  -h  .OOl(e^)*),  whidi  is  near  a  local  solution. 
Observe  how  the  diagnostics  OTl  and  OT2  tend  to  zero  at  essentially  quadratic  rate. 


5 


The  same  holds  for  the  true  error  ||  /*  — /*  ||l9-  Compare  Table  1  (Newton’s  method 
-  quadratic  rate  of  convergence)  with  Table  2  (previous  method  Disk  iteration  -  linear 
rate). 


Table  1:  levtoa  iteration  run 


It 

Value 

Optiamllty  Tasts 

Error 

SupGanaa 

*  • 

•  • 

OTl 

0T2 

■ED 

0 

3.8165288771S11E«03 

:  5.7E-01  : 

8.3E-03 

• 

1.6E400  : 

I/A 

1 

3 . 800917790397SE>03 

:  S.3E-02  : 

4.9E-04 

• 

3.9E-01  : 

8.  E-13 

2 

3 . 80008 12563274E403 

:  2.7E-03  : 

4.  E-06 

• 

1.SE-02D: 

l.lE-11 

3 

3 . 8000001798719E<K)3 

:  4.4E-06  : 

6.3E>08 

2.  E-05C: 

1.9E-11 

4 

3 . 8000000000006E>03 

:  l.SE-11  : 

2.4E-13 

3.6E-10D: 

S.3E-10 

Table  2:  Disk  iteration  nethod  run 


It  :  Value 

SupCassM 


0  :  3.816S268771SllE-»03 

1  :  3.8014192462809E+03 

2  :  3.8004689130398E<»03 

3  :  3.80011606S3124E403 

4  :  3.8000022236943E-»03 

5  :  3.8000007986885E403 

6  :  3.800000i592016E<H)3 

7  :  3.80000006292S6E«03 

8  :  3.8000000028458E403 


:ilfk-f*ll:  Optiaalitj  Tests  Error 


:  S.7E-01 
:  1.4E-01 
;  1.3E-02 
:  4.6E-93 
:  3.  E-03 
:  l.lE-03 
:  7.1E-04 
:  1.3E-04 
:  l.lE-04 


OTl 

0T2 

lED 

8.3E-03 

1.6E«00 

I/A 

5.6E>04 

1.1E<H}0D 

3.1E<H)0 

3.  E-04 

5.4E-020 

l.iE-02 

6.2E>0S 

4.  E-02C 

l.lE-03 

9.8E-07 

3.3E-020 

8.4E-05 

4.5E-07 

1.1E-02C 

3.3E-06 

1.2E-07 

7.6E>03D 

1.5E-06 

3.7E-08 

1.4E-03D 

8.7E-07 

1.9E-09 

1.2E>03C 

4.9E-0S 

Time  domain  constraints 

Recently  we  were  able  to  add  time  domain  constraints  to  (OPT)  and  obtain  optimality 
conditions  extending  Result  4  to  this  case.  We  consider  a  constrained  optimization  proUem, 
named  Constr>OPT,  where  the  minimization  is  done  over  analytic  functions  (/i./a)  that 
satisfy  a  fpven  set  of  constraints 

/  /iGi^dd  +  f  >  0,  /=:l,...,n 

Jo  Jo 

where  the  functions  are  analytic.  We  obtained, 
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Result  6  Given  F  a  smooth  function  and  constraints  as  above.  Necessary  an  sufficient  for 
a  smooth  function  T*  in  satisfying  o(e**)  =  »  never  0  on  T  to  be  a 

local  solution  to  Constr-OPT  is 

I  r(e*',T*(e‘*))  is  constant  in  e". 

II  There  exist  Fi  and  Ft  analytic  on  the  disk,  A  a  positive,  function  on  the  circle,  and 
nonnegative  constants  >  0, . . . ,  Kn  such  that  for  all  e'*  £  T, 

(e**,r*(e'*))  =  A(e**)  (e**Fi(e^)  +  +  . .  •  +  «nGi,n) 

^  (e",T*(e‘*))  =  A(e'«)  (e‘*Fa(e‘*)  +  kiU^i  +  •  •  •  + 

III  A  condition  on  second  derivatives  of  F. 

Further  analysis  shows  that  Results  4  and  6  mesh  very  well  for  the  purpose  of  construct¬ 
ing  computer  algorithms.  We  have  worked  out  such  algorithms  and  began  testing. 

Of  independent  interest  is  that  all  of  this  represents  a  new  connection  between  engineer¬ 
ing  and  an  existing  branch  of  the  mathematics  area  Several  Complex  Variables. 

Computer  algebra  for  systems  research 

There  has  been  substantial  work  on  computer  algebra  for  enj^eering  problems.  For 
example,  one  specifies  the  systems  or  circuit  components  as  letters  say  Ri  Ri  Ci  Ci  for 
resistor  and  capacitor  values  and  the  computer  produces  the  formula  for  the  transfer  function 
(no  matter  how  formidable).  Then  one  can  manipulate  it  on  the  computer. 

Our  approach  is  quite  different.  If  one  reads  a  typical  article  on  A,B,C,D  systems  in  the 
control  transactions  one  finds  that  most  of  the  algebra  involved  is  non  commutative  rather 
than  commutative.  Thus  for  symbolic  computing  to  have  much  impact  on  linear  systems 
research  one  needs  a  program  which  wiU  do  noncommuting  operations.  Mathematica,  Mac- 
syma  and  Maple  do  not  (contrary  to  what  a  salesman  will  teQ  you).  For  example,  the  most 
basic  command 

Expand[A  *  -b  C)] 

gives  A  *  B  +  A  *  C  if  A,B,C  commute  but  not  if  they  do  not.  We  have  a  package 
NC Algebra  which  runs  under  Mathematica  which  does  the  basic  operations,  block  matrix 
manipulations,  and  other  things.  The  package  might  be  seen  as  a  competitor  to  a  yellow 
pad.  Like  Mathematica  the  emphasis  is  on  interaction  with  the  program  and  flexibility. 

Mins  and  maxes  of  hamiltonians  Originally  we  wrote  the  package  to  do  linear 
control  research.  In  particular,  the  main  object  in  studying  CTRL  is  the  an  energy  balance 
(game  theoretic)  hamiltonian.  For  linear  plants 

H{x,  z,  W,  Vs)  =  V.s(x,  zf{Ax  +  BiW  +  Bicz) 

-W'^W  -I-  \\Cxx  -I-  Dx2Cz\\^  +  V,€(i,  x)^l6(Cax  -»■  D^xW)  +  ox] 
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In  the  notation  of  NC Algebra  it  is 

Ham  =  iTp{GEx[x,  z]]**(A***  +  Bl  **W  ^-82  **c  **z)-  tf[W]  **W 
(*P{*1  *  * *  * *P{®1  ♦  * f?12)  ♦  ♦  (Cl  *  ♦  *  +  D\2  **e  **z) 
+t^GEzlx,  z]]  ♦  ♦  (6  ♦  ♦  (C2  ♦  ♦  X  +  D21  *  *  W)  +  o  **z) 

One  must  compute  critical  points  (maxes  or  mins)  of  this  in  H',  a,  ft,  c  in  various  orders  which 
of  course  while  routine  is  a  tedious  process.  Also  any  variation  on  the  problem  produces 
a  new  hamiltonian  and  requires  another  tedious  computation.  NCAlgebra  automates  this. 
For  example, 


critW  =  Crtt(JJam,  W];  HnoW  =  Ham/.eritW; 
critWc  =  Crit[HnoW,^\  HnoWe  =  HnoW /  .eritWc\ 

finds  the  critical  point  of  Ham  in  W  then  in  c  and  and  evaluates  Ham  at  these  critical 
points.  This  same  4  lines  applies  to  hamiltonians  which  arise  in  other  control  problems. 

Simpliflcation  of  messy  formulas  While  NCAlgebra  can  be  used  as  a  yellow  pad  we 
are  be^nning  to  add  serious  automatic  simplification  commands.  Wavrik,  Stankus  and  I 
are  now  doing  research  in  computer  simplification  for  A,  B,  C,  D  type  linear  systems,  in  a 
highly  noncommutative  setting.  The  objective  is  in  each  particular  situation  to  find  a  list 
of  simplifying  rules.  A  complete  list  of  rules  (called  a  Grobner  basis  GB)  has  the  prc^rty 
that  if  it  is  applied  to  an  expression  until  nothing  changes  then  the  expression  is  as  simple 
as  possible  in  a  certain  sense.  Recently,  Wavrik  and  I  obtained 

Result  7  For  the  formulas  which  occur  in  studying  energy  conserving  (lossless)  systems. 
The  GB  while  infinite  can  be  summarized  as  a  list  of  Si  rules  some  of  which  depend  on  an 
integer  parameter.  We  give  the  list.  It  is  a  powerful  tool  for  studying  a  particular  class  of 
systems.  The  list  was  discovered  last  year  and  actually  proved  (with  Stankus)  to  be  a  GB 
very  recently. 

A  subset  of  these  rules  is  now  in  a  function  NCSimplifyRational{  expression]  inside  our 
NCAlgebra  package.  They  are  very  effective  on  a  limited  class  of  expressions  but  even  that 
makes  them  very  useful.  Now  we  have  some  experience  in  areas  which  use  Lyapunov  and 
Riccati  equations  and  a  line  of  experiments  involving  systems  theory  computations  which 
explore  them.  It  is  becoming  dear  that  this  is  tricky  business  and  we  are  developing  strate¬ 
gies  to  use  computer  algebra  to  obtain  systenu  theory  results.  This  is  »  matter  of  putting 
equations  in  a  simplified  canonical  form.  This  area  is  wide  open  since  in  noncommuting 
situations  the  implications  for  linear  systems  theory  are  not  explored. 

Technology  IVansfiBr 
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We  have  two  computer  programs  which  run  under  Mathematica  which  are  publically 
available. 

NC Algebra  our  non  commuting  algebra  package  has  potential  applications  in  many 
fields.  Recently  Mathematica’s  mathsource  started  distributing  it  and  in  fact  they  appear 
to  be  recommending  it  widely. 

OPTDesign  our  classical  control  program  is  available  from  ns  (send  request  to  anoptiiosiri$.ucsd.edu). 
We  do  not  intend  to  start  pushing  it  heavily  until  our  book  is  finished,  since  this  is  the  only 
account  which  ties  everything  together.  Recently,  we  completed  a  major  cleaning  of  the 
most  basic  version  of  Anopt,the  optimization  program  underlying  OPTDesign.  This  was  to 
prepare  for  porting  it  to  Matlab,  which  other  groups  have  expressed  an  interest  in  doing. 

Another  level  of  transfer  is  from  pure  to  applied  mathematics.  For  example,  in  the  last 
decade  progress  in  control  was  expedited  by  dose  connections  with  operator  theorists 
who  were  oripnally  in  pure  mathematics  but  who  now  work  the  mathematics  of  engineering 
systems.  This  originated  with  discoveries  by  DeWilde,  Fuhrmann  and  I  made  in  the  early 
1970’s. 

The  work  on  optimization  over  analytic  functions  represents  a  new  connection  between 
engineering  and  an  existing  branch  of  several  complex  variables.  Now  little  collaboration  ex¬ 
ists  between  workers  in  these  areas.  A  bi-product  of  our  development  of  (OPT)  is  possiblely 
that  a  new  group  of  pure  mathematidans  will  become  interested  in  engineering. 
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V.  Professional  Talks  &  Presentation  of  Papers 

1991-1994 


Japan -June  1991 
Workshop  on  operator  theory. 

Speaker  and  the  US  organizer. 

Sapporo  “  Hokiado  Univ. 

Japan -June  1991 

Mathematical  Theory  of  Networks  and  Systems  conference: 

Three  talks. 

Kobe 

Israel-  March  1992-sponsored  by  Binadonal  Foundation 
Conference  on  operator  theory 
Speaker 
Beer  Shev 

Israel-  March  1992-sponsored  by  Binadonal  FoundationConference 
on  operator  theory 
Speaker 
Tel  Aviv 

Southern  Calif  Partial  Differential  equations  Conference. 

April  10  -12, 1992 
Speaker 

University  of  California  Santa  Barbara 

Great  Plains  Operator  Theory  Symposium  -  May  14  -  17, 1992 
Pleniary  address 
University  of  Iowa 

Non  linear  systems  confence  -  May  29-30, 1992 
Speaker 
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Washington  University,  St  Louis  ,Mo. 

Tiawan-  sponsored  by  Tiawan  Science  Foundation  -  July  8-10,  1992 
like  a  NSF  regional  conference  in  the  US 
7  hour  lecture  set 
Tunghai  Univ. 

Conference  on  the  Interface  of  Math  and  physics  -  July  12-14, 1992 

Academica  Sinica 

Speaker 

Taipei 

Chinese  University  Hong  Kong,  July  1992 

Talk  canceled  due  to  hurricane  v>^ch  shut  down  the  dty. 

Hong  Kong,  CHINA 

SIAM  Control  Conference  On  central  organizing  committee. 

September  17  - 19, 1992 

Institute  Mathematical  Applications 

Raddison  Hotel  Minneapolis 

Minneapolis,  Minnesota 

H.  Widom's  60th  birthday  Conference ,  September  20  -  22, 1992 
University  of  California,  Santa  Cruz 
A  principal  speaker. 

Santa  Cruz,  CA 

In  residence  at  Institute  for  Mathematical  Applications, 

October  7  -16, 1992 
Special  year  in  control. 

Was  invited  to  spend  the  year  but  I  opted  for  2  weeks. 

Minneapolis,  MN 

Harvard  Mathematics  Department,  October  26-29, 1992 

Seminar 

Harvard,  MA 

MTT  October  26-29, 1992 

Seminar 

Cambridge,  MA 

Courant  Institute  of  Mathematical  Sciences,  November  2, 1992 
Colloquim 


New  York,  NY 


Virginia  Polytechnic  Mathematics  Department 
November  3  -  20, 1994 
Seminar  Math  Dept 
Blacksburg  Virginia 


Conference  on  Decision  and  Control,  December  14-19, 1992 
Presented  papers  with  J.  Ball;  with  T.  Walker  &  O.  Merino; 

&  with  W.  Zhan 
Tuscon,  Arizona 

American  Mathematical  Society  Meeting,  January  1993 
Special  Session. 

San  Antonio,  TX 

Conference  at  Free  University  of  Amsterdam  July  1993 
Speaker 

Amsterdam,  The  Netherlands 

Workshop:  Operator  Theory  and  Boundary  Eigenvalue  Problems 

July  27  -  30, 1993 

Speaker 

Vienna  Technical  University 
Wien,  Austria 

Int'l  Symposiiun:  Mathematical  Theory  of  Networks  and  Systems 
August  2-6, 1993 
Pleniary  Address 
University  of  Regensburg 

Workshop:  Algebra  and  Networks,  August  8  -15, 1993 
Speaker 

Val-Cenis,  FRANCE 

Conference  on  Decision  Control,  Decembr  1993 
Presented  papers  with  W.  Zhan;  &  with  O.  Merino 
San  Antonio,  TX 

American  Mathematical  Society,  January  1994 
Special  session  speaker 
Cincinnati,  OH 
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